Bose-Einstein condensation with a finite number of particles in a power-law trap 



A. Jaouadi, 1,2,3 M. Tclmini, 3,4 and E. Charron 1 ' 2 

1 Universite Paris-Sud, Institut des Sciences Moleculaires d'Orsay (ISMO), 91405 Orsay cedex, France 

2 CNRS, 91405 Orsay cedex, France. 
3 LSAMA, Department of Physics, Faculty of Science of Tunis, 
University of Tunis El Manar, 2092 Tunis, Tunisia 
4 National Centre for Nuclear Science and Technology, Sidi Thabet Technopark, 2020 Ariana, Tunisia 

(Dated: February 28, 2011) 

Bose-Einstein condensation (BEC) of an ideal gas is investigated, beyond the thermodynamic 
limit, for a finite number N of particles trapped in a generic three-dimensional power-law potential. 
We derive an analytical expression for the condensation temperature T c in terms of a power series 
in xo = so/fesTc, where £0 denotes the zero-point energy of the trapping potential. This expression, 
which applies in cartesian, cylindrical and spherical power-law traps, is given analytically at infinite 
order. It is also given numerically for specific potential shapes as an expansion in powers of xo up to 
the second order. We show that, for a harmonic trap, the well known first order shift of the critical 
temperature AT C /T C tx TV -1 / 3 is inaccurate when N ^ 10 5 , the next order (proportional to N^ 1 ^ 2 ) 
being significant. We also show that finite size effects on the condensation temperature cancel out in 
a cubic trapping potential, e.g. V(r) cx r A . Finally, we show that in a generic power-law potential of 
higher order, e.g. V(r) tx r a with a > 3, the shift of the critical temperature becomes positive. This 
effect provides a large increase of T c for relatively small atom numbers. For instance, an increase of 
about +40% is expected with 10 4 atoms in a V(r) tx r 12 trapping potential. 

PACS numbers: 03.75.Hh, 05.30.Jp, 64.60.-i, 37.10.Gh 



I. INTRODUCTION 

Bosc-Einstcin condensation, first predicted in the years 
1924-25 [1, 2], was observed in 1995 on rubidium and 
sodium vapors [3, 4]. These remarkable experiments have 
generated since then a clear interest in the critical prop- 
erties of this dilute macroscopic quantum state of matter. 
In particular, the transition temperature T c is a crucial 
parameter whose exact determination has been a matter 
of discussion for decades [5, 6]. 

It is indeed rather difficult to estimate the influence 
of atomic interactions on the condensation temperature 
for a uniform dilute weakly interacting Bose gas. This 
comes from the fact that this is a many-body problem 
affected by long-range critical fluctuations which have to 
be described non-perturbatively [5-7] . As a consequence, 
many different predictions were produced since the late 
1950's, like an increase or even a decrease of T c propor- 
tional to y/d7^, to a s , to a s or even to a s lna s , where 
a s denotes the atomic s-wave scattering length [6]. It is 
only in 1999 that this question was settled with the rig- 
orous demonstration of the linear increase of the critical 
temperature with a s [8]. It was later shown [9] that the 
relative correction AT C /T® behaves at second order as 
ciS + (c' 2 In S + c 2 )5 2 , where Tj? denotes the condensation 
temperature of the ideal gas in the thermodynamic limit, 
and 5 = p 1 / 3 ^, p being the (uniform) atomic density. In 
this expression, c\ ss 1.32, c' 2 w 19.75 and c 2 ~ 75.7 [7]. 

On the other hand, in a harmonic trap long-range crit- 
ical fluctuations are suppressed, and the leading (linear) 
order correction to T c can be calculated by simple pcr- 
turbative methods [10]. Using a non-perturbative ap- 
proach, it was shown in 2001 that in such a potential 



the relative correction AT C /T® behaves also at second 
order as C\5' + (c' 2 In 6' + c 2 )S' 2 , with 5' = ci s /\t, \t be- 
ing the thermal de Broglie wavelength [11]. In this case, 
ci « -3.426, c' 2 « -45.86 and c' 2 ' « -155.0 [7]. 

These results were then extended to power-law poten- 
tials by Zobay et al in a series of beautiful papers first 
using mean-field theory [12], then using renormalization 
group theory [13, 14], and later on using variational per- 
turbation theory [15]. These different theoretical studies 
helped bridging the gap separating homogeneous from in- 
homogeneous potentials by showing how an increasingly 
inhomogeneous potential suppresses the influence of crit- 
ical fluctuations. In these studies, the shape of the poten- 
tial could be summarized by a simple parameter rj which 
varies between i for the uniform gas (homogeneous case) 
and 2 for the (inhomogeneous) harmonic trap. In Rcf. 
[15], it was shown for instance that critical fluctuations 
have a marginal impact on the transition temperature 
when -q > 0.7 for At 3> a s . 

The amplitude of the corrections due to atomic inter- 
actions increases of course with the number of particles. 
On the other side, independently of atomic interactions, 
with a finite number of particles the transition tempera- 
ture is different from the one obtained in the thermody- 
namic limit. In the case of a harmonic trap, this finite 
size effect results in a downward shift of the transition 
temperature, which is, at lowest order, linear in the zero 
point energy of the trapping potential [16-18]. 

Following these developments, the present paper deals 
with higher order finite size corrections in a harmonic 
trap and in a generic power-law potential. In our ap- 
proach, the trapped Bose gas is described within the lo- 
cal density approximation. We also adopt a perturbativc 
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TABLE I. Trap parameters q and m, with i £ {1, ...,q}, and associated trapping potential U(r). The lines (a), (b) and (c 
correspond to cartesian, cylindrical and spherical power-law traps, respectively. 

q ni n 2 n 3 U(r) 

(a) cartesian trap 3 111 Ui \x/di \ ai + U 2 \y/d 2 \ a2 + U 3 \z/d 3 \ a3 

(b) cylindrical trap 2 2 1 - Ui (p/di)" 1 + U 2 \z/d 2 \ a ' 2 

(c) spherical trap 1 3 - - Ui [r/d\) ai 



approach which is only valid in the case of an inhomo- 
geneous trap. Indeed, a non-perturbative treatment is 
necessary for an accurate estimation of T c in the case of 
the uniform gas [19]. 

The plan of the paper is as follows. In section II, we 
introduce our theoretical model. We derive a general ex- 
pression for finite size effects in section III. The next 
section IV is devoted to the discussion of a specific ap- 
plication in the case of atoms confined in a crossed blue- 
detuned optical trap. In the final section V, we present 
a short summary of our main results. 



II. THEORETICAL MODEL 

A. Generic trapping potential 

For the sake of universality, and in order to treat carte- 
sian, cylindrical, and spherical traps on an equal footing, 
we follow the approach chosen in Ref. [20]. We thus 
consider a degenerate Bose gas trapped in a generic 3- 
dimcnsional power-law potential 



t/(r)=I> 



(1) 



where r^, with i £ {1, . . . , g}, are the q radial coordinates 
in the rii-dimensional subspace of the 3-dimensional space 
(see Table I for details). [/, and d, are energy and length 
scales associated with the trap. The sub-dimensions rij 
obviously verify 



3. 



(2) 



Table I gives, with the practical expression of the trap- 
ping potential U(r) of Eq. (1), a summary of the different 
parameters involved in cartesian, cylindrical and spheri- 
cal coordinate systems. 

In order to study finite size corrections to the ideal 
condensation temperature Tj? obtained in the thermody- 
namic limit, we first define the volume of the trap V(e) 
which includes the portion of space where U(r) ^ e. The 
variable e denotes here an arbitrary energy. Formally, 
this volume can be written as 



V(e) = 



d 3 7 



(3) 



For an isotropic harmonic trap of angular frequency w, 
this definition gives a volume of at the energy e = 

iftw, where et^ is the usual characteristic length of the 
harmonic trap. In the general case, integration over the 
3-dimcnsional space yields the expression 



V(e) 



en 
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where 
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(4) 



(5) 



and where 9 = 8, 2vr or 47r/3 in the case of a cartesian, 
cylindrical or spherical trap respectively. In these equa- 
tions, T(x) denotes the complete gamma function [21]. 
In addition, in Eq. (4), r\ denotes a very important pa- 
rameter which characterizes the shape of the potential. 
This parameter is defined by 



E 

i=l 



(6) 



One can see in this last equation that the shape pa- 
rameter r\ equals 2 in a harmonic trap corresponding to 
on = 2 (Vz), while 77 equals ^ for a square- well poten- 
tial corresponding to a, -> 00 (Vi). Intermediate values 
h < rj < 2 correspond to the different cases of power-law 



potentials, as defined in Eq. (1). 



B. Density of states 

In the theoretical description of Bose-Einstein con- 
densation within the local density approximation, three- 
dimensional integrals of functions of the trapping poten- 
tial U(r) are common. They can be greatly simplified 
using the following transformation [12] 



f(U(r))d 3 r = / f(e)v(e)de, 



(7) 



where the spatial volume v(e)de of the equi-potential 
shell of U (r) with width de at energy e can be obtained 
from Eq. (4). This yields 



v(s) 



en 



(8) 
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The transformation (7) is for instance helpful for the de- 
termination of the density of states 17(e) associated with 
the generic trapping potential U(r). In a semi-classical 
approach, this density is given by 



9(e) 



1 



d 3 r /// d 3 p<5(e-e cl (r,p)), (9) 



where p is the momentum vector of the atom and where, 
in the case of the ideal gas, e c i(r,p) denotes the disper- 
sion relation 
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Radial coordinate 



£cl (r, P ) = ^ + C/(r) 



(10) 



m being the atomic mass. In the expression (9), S(x) 
denotes the Dirac delta function. Since e c i(r,p) does 
not depend on the orientation of the momentum vector 
but on its magnitude only, one can integrate over the 
momentum coordinate to find [221 



l(e) 
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[e-U(r)] 2 d 3 r. 



(11) 



Integrating this expression using the transformation (7) 
then yields 



_ 1 /m\! en 



(12) 



FIG. 1. (Color online) Schematic representation of a one- 
dimensional harmonic trapping potential V(x) oc x 2 (solid 
blue line in the left panel) and of a one-dimensional power- 
low trapping potential V(x) oc x 12 (solid blue line in the right 
panel), with their associated energy level structures. 



where the spatial density distribution of the thermal gas 
is given, in the local density approximation, by 



1 



nr(r) = T3- 9# ( exp 



fi - U(r) 



knT 



(14) 



In this expression /j is the chemical potential and Xt the 
thermal de Broglic wavelength at the absolute tempera- 
ture T 



Consequently, g(e) is proportional to e 2 for a harmonic 
trap and to yfe in the homogeneous limit corresponding 
to 77 —> i . This strong evolution of the density of states as 
a function of the shape of the potential could already be 
inferred from Fig. 1, which shows a schematic represen- 
tation of a one dimensional harmonic potential (left side) 
and of a one-dimensional power-law potential V(x) oc x 12 
(right side), with their associated energy level structures. 
These two traps of different shapes are characterized by 
very different energy level structures, and are therefore 
characterized by very different densities of states. Since 
the condensation process consists in an amplified growth 
of the ground state population resulting from stimulated 
collisions of atoms in a thermal reservoir, these very dif- 
ferent densities of states lead to different condensation 
temperatures [23], but also, as we will see hereafter, to 
different corrections due to finite size effects. 



At = 



/ 27Tft 2 



(15) 



g s (z) denotes the Bose function of order s, which, for 
our purposes, can be considered as defined in terms of its 
series representation 



9s(z) 



Y- 

fc=l 



(16) 



In the case of an ideal gas in the thermodynamic limit, 
the critical value of the chemical potential is \i = 0. In- 
troducing this value in Eq. (14) and integrating Eq. (13) 
using the transformation (7) yields 



C. Ideal condensation temperature in the 
thermodynamic limit 

Assuming that, at the condensation temperature, the 
number of condensed atoms is negligible, one can evalu- 
ate the total number of particles N as 



N = 



n T (r)d 3 r, 



(13) 



N 



en 

If 
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ri+l 



(17) 



from which the critical temperature T° of the ideal gas in 
the thermodynamic limit can be obtained. T° therefore 
varies as N 1 ^ 71 ^, the shape parameter 77 being defined 
in Eq. (6). Since \ < 77 ^ 2, we see already that the 
shape of the potential has a strong impact on the depen- 
dence of the critical temperature Tj? with the number of 
particles [23]. 
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III. FINITE ATOM NUMBER CORRECTIONS 



and for r\ <E N, one obtains 



A. Derivation of the shift 

For a finite number of ideal particles, the zero-point 
energy of the trapping potential can not be neglected 
anymore and the critical value of the chemical potential 
is bounded by fj, = eq. Following the method used pre- 
viously to derive Eq. (17), one can introduce this value 
of /j in Eq. (14) and integrate Eq. (13) using the trans- 
formation (7) to obtain a new expression for the total 
number of atoms which is very similar to Eq. (17): 



N = 



en 



/ m \ ■ 
\2tt) 



(v\xo)l 



J7+1 



(18) 



where 



where 



T(r, 



''() ! 



ln(x ) 



^ k 

k=l 



(26) 



(27) 



the parameters Bk, v and Ck, n being already defined in 
Eqs. (25b) and (25c). 

Practically, for small values of xq, one can truncate the 
sums in Eqs. (24) and (26) to the desired level of accuracy 
and introduce the expression obtained for (r]\xo) x in Eq. 
(20) to get the correction AT C = T c - T c ° due to the 
presence of a finite number of particles. 



(v\xo) a 



gsie 



x v 2 dx , 



(19) 



and where xq = So/ksTc. Note that the change of vari- 
able x = e/ksTc was used. The condensation tempera- 
ture T c = T c ° + AT C now includes the corrections induced 
by the presence of a finite number of particles. Compar- 
ing Eqs. (17) and (18), we find 



J7+1 



(v\xo)Z 



(20) 



The calculation of the integral (^|^o)^ requires to sep- 
arate long from short range effects using 



(21) 



The first integral is easily obtained using the Bose func- 
tion (16) as 



while the second integral can be expressed as 



(v\xoT ° 
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(>o(l- 



du . 



(22) 



(23) 



In these last two expressions, one can expand the two 
Bose functions g v +i(e Xo ) and 173 (e 2 ^ 1 " 11 )) using the se- 
ries representations of these functions about zero [24] . 
For non- integer values of r\ (i] G R\N), this yields 



(v\xo)l 



A 



OO 



Bk, v — C, 



k,r/ X Q 



k=0 



b • 



where 



Bk,rj 
Ck,rj 



i)cos(?77r)r(-?7) 
i)C(r? + l -*)/*! 

i)C(i-fc)/rfa + * + *), 



(24) 



(25a) 
(25b) 
(25c) 



B. Discussion 

1. The harmonic trap 

In the special case of a harmonic trap (77 = 2), and for 
xq ^0.1, the truncation can be limited to second order in 
xq, and a simple expansion of T c /T® [Eq. (20)] in powers 
of xq yields at second order 



AT C 

x c 

where 



C(2) 



;Xq 



C(l) 



3C(3) 3C(3)r(f) 



Xq - a(x a )xl , (28) 



a(x ) 



3C(3) 



2[C(2)]- 
3C(3) 



3 
4 



ln(xo) 



(29) 



In this expression, one can recognize the well-known first 
order correction — (£(2)/[3£(3)]) xq ~ — 0.4561 xq of the 
harmonic trap [25] . Note however that the next order cor- 
rection + (C(|)/[3C(3)r(|)])4 /2 ~ +0.5449 x 3 /2 should 

be taken into account when xo 10 -3 for an accurate 

3/2 

estimation of the correction to T c . Up to the order x , 
the correction (28) to the condensation temperature in 
a harmonic trap can be rewritten as a function of the 
number of atoms as 



where 



and 



C(2)cV*-* +C(| / 2 ^ 



37^(3)^ 



7Y"2 



2C(3)*a; 

-0.7275 (^) N~i + 1.098 f 



= - THUi 



n 



(31a) 



(31b) 
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TABLE II. Practical second order expressions of the shift to 
T c due to finite size effects for different values of the shape 
parameter 77. T° denotes the condensation temperature ob- 
tained in the thermodynamic limit and xo = eo/ksTc denotes 
the ratio of the zero point energy of the trapping potential to 
the thermal energy ksTc. 
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AT C /T° 



3/4 
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+0.3037a;o 



0.8393 Xo' " + 0.9685 x\ /2 + 0.2239 x° 



3/4 



0.2921 x 7 /4 + 0.3222 xf, + 0[xl 



0.2536 4 /4 + 0.2064 xl' 2 



9/4, 



FIG. 2. (Color online) Relative correction AT C /T C ° due to 
finite size corrections as a function of the number of atoms 
N (logarithmic scale) in an isotropic harmonic trap. The 
corrections up to order iV~ 1//3 and AT -1 ' 2 in Eq. (30) are 
shown as a solid line with black circles and a dashed line with 
red squares, respectively. 



denote the arithmetic and geometric means of the oscilla- 
tor frequencies oji associated with each radial coordinate. 

This relative correction AT C /T® in shown in Fig. 2 
as a function of the number of particles N in the case 
of an isotropic harmonic trap (w m = uj). In this figure, 
the solid black line (with circles) represents the relative 
correction obtained using the first order correction only, 
whereas the dashed red line (with squares) represents 
the relative correction obtained using the more accurate 
expression (30). 

One can see in Fig. 2 that the well known first order 
correction — 0.7275 (w m /ZU) N~ 3 is inaccurate for small 
atom numbers. This is due to the fact that, beyond the 
first order correction, the next order is not a second or- 
der correction, which would be proportional to N~ 3 , but 
rather a correction proportional to N~~ 2 . When using 
only the first order correction found in the literature, the 
error made on the shift of the critical temperature is of 
about 91%, 48%, 29% and 18% for N = 10 3 , 10 4 , 10 5 
and 10 6 , respectively. We can therefore conclude that 
the last term of Eq. (30), proportional to iV~2, should 
be included for an accurate estimation of finite size cor- 
rections to T c in a harmonic trap when N ^ 10 5 . 



2. The power-law trap 

Beyond the simple case of the harmonic trap, a general 
expansion of T c /T® [Eq. (20)] in powers of xo can be 
derived for the two cases 77 € R\N and 77 € N. We give in 
Table II, at second order, practical expressions for a set 
of specific values of 77 between 0.75 and 2, corresponding 
to different kinds of inhomogcneous potentials. 

One can first notice in Table II that the case 77 = 3/2 is 
a special case. Indeed, for 77 = 3/2 all finite atom number 
corrections cancel out in Eq. (24) except for the constant 



0.8960 x\ /2 + [0.9003 + 0.3040 In x ] x Q 
1 +[0.6474 + 0.8171 In x ]x 3 /2 

+0.139(0.201 + lnxo] [11.01 + lnsco] xl + O[x 5 /2 ] 



5/4 


0.8652 x 3 /4 - 1.399 x + 0.8440 x h J 4 
+ 1.216 xl /2 - 4.209 x 7 J 4 + 5.675 xl + O[x ] 9/4 


3/2 


(to all orders in xo, see text for details) 


7/4 


-0.5662 x + 0.6653 xf/ 4 - 0.5636 x 7 Q /4 
+ 1.098 xl + 0[xl /4 ] 


2 


-0.4561x0 + 0.5449 x 3 J 2 
+[0.2082 + 0.1387 lnx ] ^0 + 0[xf/ 2 ] 



term B s = C(|)- T c is then identical to T r ° for all N. 
This cancellation comes from the fact that incidentally 
Cfe,3/2 = -Bfc+1,3/2 and A 3/2 = [see Eqs. (25)]. This is 
the only value of the shape parameter 77 where there is rig- 
orously no finite atom number correction, whatever the 
number of trapped atoms. This special case corresponds 
for instance to an isotropic cubic trap, with V(r) oc r , 
but also to a trapping potential of cylindrical symmetry 
such as V(r) oc p z + |z| 3 or to V(r) oc |x| 3 + |y| 3 + |z| 3 
in the case of a cartesian power-law trap. It also applies 
for instance to V(r) oc p A + z 2 , or to any combination of 
power- law potentials for which 77 = 3/2 [see Eq. (6)]. 

The variation of the relative correction AT C /T° is 
shown MS i\ function of the shape parameter 77 in Fig. 3 
for xo = 0.1 (solid line with black circles) and xo = 0.01 
(dashed line with red squares). From the magnified in- 
set of this figure, one can notice again that the point 
77 = 3/2 plays a particular role since trapping potentials 
with 77 > 3/2 are characterized by a decrease of T c due 
to finite size effects, while 77 < 3/2 leads, on the contrary, 
to an increase of T c . This could already be inferred from 
Table II, where it is seen that the lowest order correction 
is positive when 77 < 3/2 and negative when 77 > 3/2. 

In addition, one can notice in Fig. 3 that, even when 
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FIG. 3. (Color online) Relative correction AT C /T° (in %) due 
to finite size effects as a function of the shape parameter r\ 
for xq = 1CP 1 (solid line with black circles) and xo = 1CP 2 
(dashed line with red squares). The inset is a simple mag- 
nification of the region surrounding r\ = 3/2. Note that 
AT C /T C ° = for n = 3/2, Vz . 



xq <C ksT c , the increase in T c becomes significantly large 
when 77 ^ 1. Indeed, for r) = 3/4, an increase of 37% 
is expected when xo = 0.01, and an increase of 84% is 
expected when xq — 0.1. A large upward shift of T c 
can therefore be expected in the case of a generic high- 
order power-law potential, like V(r) cx r a with a > 3 for 
instance. One should note here that we have limited our 
study to 0.75 ^ 77 sC 2, i.e. to inhomogeneous potentials 
since a non-perturbative treatment would be necessary 
for an accurate estimation of AT C in the homogeneous 
limit 77 — > \. 

In order to check the influence of these corrections in 
a realistic trapping situation, we now turn to the case 
of a rubidium gas trapped by blue-detuned Laguerre- 
Gaussian optical beams. 



IV. APPLICATION TO BLUE-DETUNED 
LAGUERRE-GAUSSIAN OPTICAL TRAPS 

A. Introduction and model system 

Recently, we have proposed an original route to achieve 
Bose-Einstein condensation using dark power-law laser 
traps [23]. Experimentally, it has been demonstrated 
that cold atoms can be trapped and guided in such opti- 
cal setups [26] . These traps are created with two crossing 
blue-detuned Laguerre-Gaussian optical beams. A sim- 
ple control of their azimuthal order I allows for the explo- 
ration of a wide range of situations, with atoms trapped 
in confining power-law potentials in one, two, and three 
dimensions. It was shown that, for a fixed atom number, 
higher transition temperatures are obtained in configura- 
tions where r\ < 2, compared to a harmonic trap (rj = 2) 
of the same size. We investigate here finite size effects in 
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FIG. 4. (Color online) Relative correction AT C /T° (in %) due 
to finite size effects as a function of the shape parameter rj 
(lower scale) and of the azimuthal index I (upper scale) for two 
Laguerre-Gaussian beams of power P = 5W and detuning 
5=10 THz. The solid black line with circles corresponds to 
N — 10 3 particles, the red dashed line with squares to iV = 
10 4 and the blue dash-dotted line with triangles to TV = 10 6 . 
The volume of the trap has been chosen such that, for each 
atom number, the condensation temperature T° obtained in 
the thermodynamic limit for rj = 2 is 100 nK (see [23] for 
details). 



such optical traps. 

To create a three-dimensional trap, we consider an all- 
optical configuration consisting of two perpendicularly 
crossing Lagucrre-Gauss laser modes of same radial in- 
dex p = and azimuthal index I. The polarizations of the 
two beams are chosen to be orthogonal to avoid any in- 
terference pattern. The first beam, circularly symmetric, 
propagates along the z direction and provides trapping in 
the {x, y) plane. Trapping in the third dimension is pro- 
vided by another strongly elongated Laguerre-Gaussian 
beam shaped elliptically in the form of a light sheet and 
propagating in the x direction. If the two beams are char- 
acterized by the same laser power, waist and detuning, 
the corresponding potential near the trap center can be 
described by 



V(p,z) = U (p 2e + z 2 C) 



(32) 



Uo is then a simple function of the laser power, detuning 
and waist of the two beams [23], and 77 = (£ + 3)/(2£). 

B. Finite size effects 

The relative correction AT C /T° due to finite size effects 
in such optical traps is shown in Fig. 4 as a function of 
the shape parameter 77 (lower scale) or equivalently of the 
azimuthal index i (upper scale), for realistic values of all 
optical parameters (see figure caption). One can notice 
in this figure that, even for large condensate occupation 
numbers, the shift of the condensation temperature due 
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to finite size effects is not negligible when 77 < 1. Indeed, 
for 77 = 3/4 this correction is of about 60%, 38% and 16% 
for N = 10 3 , 10 4 and 10 6 , respectively. This potential 
corresponds to a 12 th -power power- law trap, which can 
be obtained practically using t = 6 Laguerre-Gaussian 
beams [26]. 

Azimuthal index 



• N=10 3 




0.75 1 1.25 1.5 1.75 2 



Shape parameter 

FIG. 5. (Color online) Condensation temperature in nK as 
a function of the shape parameter r\ (lower scale) and of the 
azimuthal index i (upper scale) for two Laguerre-Gaussian 
beams of power P — 5W and detuning S — lOTHz. The 
solid green line with black crosses corresponds to the conden- 
sation temperature T? obtained in the thermodynamic limit. 
The dashed black line with circles gives the condensation tem- 
perature T c for N — 10 3 particles, the red dashed line with 
squares for N = 10 4 and the blue dashed line with triangles 
for N = 10 6 . The volume of the trap has been chosen such 
that, for each atom number, the condensation temperature 
obtained in the thermodynamic limit for r\ = 2 is 100 nK (see 
[23] for details). 

Similarly, Fig. 5 shows the variation of the conden- 
sation temperature with the shape parameter 77 and the 
azimuthal index £. The temperature T° obtained in the 
thermodynamic limit is shown as a thick green solid line. 
It is identical for all atom numbers TV since we have cho- 
sen the volume of the trap such that the condensation 
temperature obtained in the thermodynamic limit is the 
same for all atom numbers (100 nK with 77 = 2, see [23] 
for details). As already discussed in [23], the temperature 
Tj? increases with £ since 

In the same figure, the condensation temperature T c tak- 
ing into account finite size effects is shown as a dashed 
black line with circles for N = 10 3 , a red dashed line with 
squares for TV = 10 4 and a blue dashed line with trian- 
gles for N = 10 6 . Compared to a harmonic trap of the 
same size, it can be concluded that much higher critical 
temperatures can be obtained in a power-law trap, and 
this especially for small condensation occupation num- 
bers. This large increase results from the combination 
of two positive effects: first the effect of the potential 
shape on the temperature T® [Eq. (17)] and second the 
finite size effect [Eq. (20)]. Indeed, with the optical pa- 



rameters chosen in Fig. 5, compared to T c ° = 100 nK in 
the harmonic case (77 = 2), the condensation tempera- 
ture for 77 = 3/4 is about 283 nK, 251 nK and 212 nK for 
N = 10 3 , 10 4 and 10 6 , respectively. This corresponds to 
an increase of T c by about 183%, 151% and 112%, re- 
spectively, when compared to the transition temperature 
found in a harmonic trap. 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have presented a detailed theoretical 
analysis of finite size effects on the transition tempera- 
ture T c of an ideal Bose gas trapped in an inhomogeneous 
generic three-dimensional power-law potential. Using the 
local density approximation, we have derived an analyti- 
cal expression of T c , expressed in terms of a power series 
in xq = e^/kBTc, where Eq denotes the zero-point energy 
of the trapping potential. We have also given numerical 
expressions truncated at second order for a specific set of 
power-law potentials. 

More precisely, we have shown that in the case of a 
harmonic trap the usual lowest-order correction is inac- 
curate when the number particles N is less than 10 5 . In 
this case, the next order has also to be taken into account, 
yielding to the relative correction 

^ ~ -0.7275 (^) J\T* + 1.098 f^f AT* , (33) 
Tg \ u J V ui / 

where ui m and ZJ denote the arithmetic and geometric 
means of the harmonic oscillator frequencies. 

Beyond the specific case of the harmonic potential, 
we have also shown that, if the atoms are trapped in a 
power-law potential characterized by a shape parameter 
77 = 3/2, all finite size effects cancel out, and the criti- 
cal temperature is simply equal to the one obtained in 
the thermodynamic limit, independently of the number 
of particles. In the case of a trap of cylindrical symmetry, 
this happens for example for V(r) oc p 3 + \z\ 3 , but also 
for V(r) oc p A + \z\ 2 . 

In addition, we have shown that in a power-law po- 
tential characterized by a shape parameter 77 < 3/2, the 
temperature shift due to finite size effects becomes pos- 
itive. This shift may be relatively large for 77 < 1 with 
relatively small atom numbers (N ^ 10 4 ). The mag- 
nitude of this phenomenon was finally confirmed in the 
specific case of rubidium atoms confined by blue-detuned 
Laguerre-Gauss optical laser beams. 

Finally, one should note that even though our approach 
gives the exact result for all types of power- law potentials 
within the limits imposed by the approximations used in 
our theoretical model, it is not able to capture the essence 
of critical phenomena taking place in a purely homoge- 
neous system, due to the importance of long-range non- 
perturbative critical fluctuations in such systems. 
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